i/1 


AD-A170  022  AN  OC/V/+/E/)  ALGORITHM  FOR  FINDING  AN  EDGE-MAXIMAL 
SUBGRAPH  UITH  A  TR-FO  (U)  CARNEG IE-MELLON  UNIV 
PITTSBURGH  PA  MANAGEMENT  SCIENCES  RESEAR  E  BALAS 
UNCLASSIFIED  JUL  85  MSRR-519  N00014-82-K-0329  F/G  12/i 


NL 


AN  0 ( I V |  +  ! E i )  ALGORITHM 


FOR  FINDING  AN  EDGE-MAXIMAL  SU3GRAPH 
WITH  A  TR-FOKMATIVE  COLORING 
bv 

Egon  Balas 


Carnegie -Mel  Ion  University 

PITTSBURGH,  PENNSYLVANIA  15213 


GRADUATE  SCHOOL  OF  INDUSTRIAL  ADMINISTRATION 

WILLIAM  LARIMER  MELLON,  FOUNDER 


- - r 


I-.r- .?  d b*  *  • 

f  b‘»  ■  •  ■  s»  n  • 

r;  \  hn_:K‘ -i. 


W. P. #8-85-86 


Management  Science  Research  Report  No.  MSRR-519 


AN  0 ( | V 1  +  |E ] )  ALGORITHM 
FOR  FINDING  AN  EDGE-MAXIMAL  SUBGRAPH 
WITH  A  TR-FORMATIVE  COLORING 
by 

Egon  Balas 


July  1985 


The  research  underlying  this  report  was  supported  by  Grant  ECS-8503192 
of  the  National  Science  Foundation  and  Contract  N00bl4-85-K-0198 
with  the  U.S.  Office  of  Naval  Research.  Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of  the  U.S.  Government. 

Management  Science  Research  Group 
Graduate  School  of  Industrial  Administration 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  15213 


ABSTRACT 


If  TR  is  the  class  of  triangulated  graphs,  a  TR-formative  edge  coloring 
is  a  green/ red  coloring  of  the  edges  of  a  graph,  such  that  the  green  graph  is 
triangulated  (i.e.  belongs  to  TR)  and  the  red  graph  has  no  triangles.  Recently 
Balas,  Chvatal  and  Nesetril  gave  an  0(|v|^)  algorithm  for  finding  a  maximum- 
weight  clique  in  any  graph  G  =  (V,E)  with  a  known  TR-formative  edge  coloring. 

In  this  paper  we  give  an  0(|v|  +  j E | )  time  algorithm  for  finding  in  an  arbitrary 
graph  an  edge-maximal  subgraph  with  a  TR-formative  coloring.  This  can  be  used 
to  construct  improved  implicit  enumeration  procedures  for  finding  a  maximum- 
weight  clique  in  an  arbitrary  graph. 

Our  algorithm  consists  of  two  subroutines,  also  of  interest  in  their 
own  right:  one  finds  an  edge-maximal  triangulated  subgraph,  the  other  one  an 
edge-maximal  triangle-free  subgraph,  in  an  arbitrary  graph. 
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1.  Introduction 

Let  TR  denote  the  class  of  triangulated  graphs,  i.e.  graphs  that  have 
no  hole  (chordless  cycle  of  length  >  4).  For  an  arbitary  graph  G  =  (V,E), 
we  define  a  TR-formative  edge  coloring  of  G  as  a  green/ red  coloring  of  the 
edges  of  G,  such  that  the  green  graph  is  triangulated  and  red  graph  has  no 
triangles.  In  a  recent  paper,  Balas,  Chvatal  and  Nesetril  [1985]  have 
given  an  0(jv|^)  algorithm  for  solving  the  maximum  weight  clique  problem 
(MWCP  for  short)  on  any  graph  G  with  a  known  TR-formative  edge  coloring. 

This  result  can  be  used  to  derive  improved  implicit  enumeration  algorithms 
for  solving  the  MWCP  on  an  arbitary  graph  G,  provided  one  has  an  efficient 
way  of  generating  subgraphs  of  G  with  a  TR-formative  edge  coloring.  Indeed, 
suppose  that  for  some  E^E,  the  subgraph  G[E7]  =  (V,  E 7)  of  G  generated 
by  E7  has  a  known  TR-formative  edge  coloring.  Then  the  MWCP  on  G[E7]  is 
solvable  in  0(|v|^)  time.  Let  K  be  a  maximum  weight  clique  of  G[E7]  with 
weight  w(K).  If  G  has  a  clique  K7  with  w(K7)  >  w(K) ,  then  G(K7)  has  at 
least  one  edge  in  E  \  E7.  Thus  an  implicit  enumeration  algorithm  can  be 
constructed  that,  as  a  branching  rule,  recursively  replaces  the  current 
graph  G  by  the  collection  of  subgraphs  G(N(e^)),  GCNCe^))  -  e^,  ..., 
G(N(e^))  -  {e^, .  .  .  ,e^_^},  where  [e^,...,e  }  =  E  \  E,7  and  for  any  edge 
e  =  (u,v), 

N(e):  =  cV \{u,v] j (w,u) tE  and  (w,v)eE), 

while  G(N(e))  denotes  the  subgraph  of  G  induced  by  N(e).  The  procedure  for 
finding  a  subgraph  with  a  TR-formative  edge  coloring  can  then  be  applied  to 
each  of  the  above  graphs,  while  any  graph  G(N(e^))  -  (e  ,...,e^  that  can 
be  shown  to  have  no  clique  of  weight  larger  than  w(K)  -  w(u^)  -  w(v  )  (where 
(u^,  v^)  =  e^),  can  be  discarded.  Naturally,  the  larger  the  subgraph  with 
a  TR-formative  coloring  that  one  is  able  to  generate,  the  fewer  branches  are 


needed;  hence  one  is  moved  to  search  for  edge -maximal  subgraphs,  where  we 
define  the  latter  as  follows.  A  triangulated  subgraph  G[F]  is  called  edge- 
maximal  (with  respect  to  set  inclusion)  if  there  exists  no  F  ^  F,  Ft  E, 
such  that  G[F;]  is  trangulated.  An  edge-maximal  triangle-free  subgraph  is 
defined  analogously.  Now  if  G[E']  is  a  subgraph  of  G  with  a  TR-formative 
edge  coloring  [F,D] ,  where  G[F]  is  triangulated,  G[D]  is  triangle-free,  and 
FUD  =  E7,  we  say  that  GfE']  is  edge-maximal  if  there  exists  no  F  F,  F  t  e\d, 
such  that  G[F7]  is  triangulated,  and  no  D1^  D,  D7  C  E\F,  such  that  G[D7]  is 
triangle-free. 

An  implicit  enumeration  procedure  of  the  above  type  was  used  by  Balas  and 
Yu  [1984]  to  find  a  maximum  (unweighted)  clique  in  an  arbitrary  graph.  The  Balas - 
Yu  algorithm  generates  a  maximal  induced  subgraph  whose  chromatic  number  is  equal 
to  the  size  of  its  maximum  clique,  and  in  which  a  maximum  clique  can  be  found 
in  0(  |v|+'(E| )  time.  If  G(S)  denotes  the  maximal  induced  subgraph  generated, 
then  G  is  replaced  by  G(N(v^)),  G(N(v2)"{v^}) , . . . ,G(N(v  )-(v^, . . . ,vp_^}) .  The 
computational  results  obtained  on  randomly  generated  graphs  with  up  to  400  ver¬ 
tices  and  30,000  edges  indicate  that  the  procedure  is  clearly  superior  to  earlier 
algorithms  that  use  straight  implicit  enumeration. 

In  this  paper  we  give  an  0(jv|+|E|)  procedure  for  finding  an  edge-maximal 
subgraph  with  a  TR-formative  coloring  in  an  arbitrary  graph.  The  procedure  con¬ 
sists  of  two  independent  algorithms,  each  of  which  is  also  of  interest  in  its 
own  right.  The  first  one  finds  an  edge-maximal  triangulated  subgraph,  the 
second  one  finds  an  edge-maximal  triangle- free  subgraph,  in  an  arbitrary  graph. 
Applied  to  a  graph  G  =  (V,E),  algorithm  I  finds  an  edge-maximal  triangulated 
subgraph  G[F]  of  G.  Applied  to  G[ENE] ,  algorithm  II  then  finds  an  edge-maximal 
triangle-free  subgraph  G[D]  of  G[E'\F].  The  resulting  graph  G[FUD]  is  then  an 
edge-maximal  subgraph  of  G  with  the  TR-formative  edge  coloring  [F,D]. 
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2.  Edge-Maximal  Triangulated  Subgraphs 
Our  algorithm  for  finding  an  edge-maximal  (EM)  triangulated  subgraph 
of  a  graph  is  based  on  the  Balas-Yu  algorithm  for  finding  a  maximal  trian¬ 
gulated  induced  subgraph.  Like  the  latter,  it  uses  the  ideas  of  a  procedure 
by  Rose,  Tarjan  and  Lueker  [1976]  for  testing  triangularity,  and  it  runs  in 
0(|v|+|e|)  time.  An  earlier  algorithm  for  finding  an  EM  triangulated  sub¬ 
graph,  by  Dearing,  Shier  and  Warner  [1984[,  requires  0(A  "1E|)  steps,  where 
A  is  the  maximum  degree. 

To  explain  the  algorithm,  we  have  to  recall  a  few  properties  of  trian¬ 
gulated  graphs.  A  vertex  is  called  simplicial  if  all  its  neighbors  are 
adjacent  to  each  other.  Every  triangulated  graph  has  a  simplicial  vertex 
(Dirac  [1961]);  it  follows  that  a  triangulated  graph  has  at  most  as  many 
cliques  as  vertices.  An  ordering  vi>*‘,>vn  the  n  =  |v|  vertices  of  a 
graph  G  =  (V,E)  is  called  perfect  if  for  i=l,...,n,  v^  is  simplicial  in 
G(fvi,v^+^, . . . »vn}) .  A  graph  is  triangulated  if  and  only  if  there  exists 
a  perfect  ordering  of  its  vertices.  Based  on  this  property,  triangulated 
graphs  can  be  recognized,  and  their  cliques  can  be  listed,  in  0(|v|  +  |e|) 
time.  (Rose,  Tarjan  and  Lueker  [1976]). 

Given  a  graph  G  =  (V,E),  our  algorithm  generates  an  EM  triangulated 
subgraph  G[F]  of  G,  along  with  a  perfect  ordering  cr  of  the  vertices  of 
G[F] .  The  ordering  a  is  generated  backwards  (i.e.  the  last  rank  is  as¬ 
signed  first,  the  next  to  last  seconds,  etc.),  based  on  a  lexicographic 

labeling  of  the  vertices.  For  some  ordering  a  =  (v. ,...,v  )  of  the  ver- 

i  n 

tices  of  G,  we  say  that  v.  is  a  successor  of  v.  if  v.  and  v.  are  adjacent 

J  -  x  j  1 

and  j  >  i.  If,  in  addition,  k  >  j  for  all  successors  v^  of  v^  other  than 

v.,  then  we  say  that  v.  is  the  first  successor  of  v, . 

J  j  -  i 

We  will  use  two  known  facts,  established  in  Rose,  Tarjan  and  Lueker 
[1976] : 


Fact  !•  If  for  i  =l,...,n,  the  first  successor  of  v  is  adjacent  in 

G[F]  to  every  other  successor  of  v^,  then  for  i  =  1 . n,  v ^  is  simplicial 

in  the  subgraph  of  G[F]  induced  by  {v. ,...,v  }. 

i  n 

Fact  2.  G[F]  is  an  EM  triangulated  subgraph  of  G  is  and  only  if 
there  exists  no  eeE\F  such  that  G[FU{e}]  is  triangulated. 

Algorithm  I  (for  finding  an  EM  triangulated  subgraph) 

0.  Initialization.  Assign  the  label  0  to  every  vertex. 

Set  i  <-  n  =  |v| ,  F  -  0,  a  -  0. 


1.  Choosing  a  vertex.  If  i  =  0,  stop:  G[F]  is  an  EM  triangulated 
subgraph  of  G  and  a  is  a  perfect  ordering  in  G[F] . 

Otherwise,  choose  an  unnumbered  (i.e.  unranked)  vertex  v  with  lexi¬ 
cographically  largest  label.  Assign  number  i  to  v,  i.e.  set  v.  *-  v,  insert 


v^  into  a  as  the  first  element,  i.e.  set 


a  -  <V  »w .  V, 


and  go  to  2. 

2.  Adding  edges  to  F.  Let  SCT(vi>  =  (w^...^)  be  the  list  of  successors 
of  v^ .  Add  to  F  the  edge  (v^.w^)  and  all  the  edges  (v^w^)  such  that  w^  is 
a  successor  of  v^ ,  adjacent  to  w^  in  G [F] ;  i.e.,  set 

F  "  F  U  (v^ , w^ )  U  {(v^,Wj)  |wj6S^(v^)  and  (w^w^jeF} 
and  go  to  3. 


3.  Labeling .  Append  i  to  the  label  of  each  unnumbered  vertex  w  adjacent 
to  v^  in  G,  set  i  •-  i  -  1,  and  go  to  1. 

Theorem  1.  Algorithm  I  generates  an  EM  triangulated  subgraph  of  G  in 
0 ( | V 1  +  |E|)  time. 

Proof.  Let  a  =  (v^,...,vn)  be  the  ordering  generated  by  Algorithm  I. 

From  Fact  1,  for  i  =  l,...,n,  v^  is  simplicial  in  the  subgraph  of  G[F]  in¬ 


duced  by  [v. ,  v  .,..., v  }.  Hence  G[F]  is  triangulated.  From  Fact  2,  to 


prove  that  G [F]  is  edge-maximal  it  is  sufficient  to  show  that  for  all 
ee£  \  F,  G[FU{e}]  is  not  triangulated.  But  this  follows  from  the  rule 
which  governs  the  addition  of  new  edges  to  F  (step  2).  Thus  G[F]  is  an 
EM  triangulated  subgraph  of  G. 

Next  we  establish  the  complexity  of  the  algorithm.  Each  application 

of  Step  3  (labeling)  requires  0(deg  v  )  operations,  and  the  list  of  labeled 

vertices  can  be  kept  in  lexicographically  decreasing  order.  Thus  Step  1 

(choosing  an  unnumbered  vertex  with  lexicographically  largest  label)  can  be 

carried  out  in  constant  time.  Step  2  (adding  edges  to  F)  again  requires 

0(deg  v  )  operations.  To  get  the  complexity  of  the  entire  algorithm,  we 

n 

sum  over  all  v^  i  =  1,  ....  n,  to  obtain  0(  £  deg  (v^)  =  o|v|  +  |e|).:; 

3.  Edge-Maximal  Triangle-Free  Subgraphs 
A  straightforward  method  for  finding  an  EM  trinagle-free  subgraph 
G [ D]  of  a  given  graph  G  =  (V,E)  is  to  initialize  the  edge  set  D  as  empty, 
then  examine  every  edge  of  E  in  some  arbitrary  order  and  put  it  into  D  if 
and  only  if  this  does  not  create  a  triangle  in  D.  Examining  every  edge  from 
this  point  of  view  requires  0 ( ) V j )  steps,  hence  the  whole  procedure  requires 
0 ( | V  |  •  |E| )  steps. 

However,  one  can  do  better  than  this  naive  approach  by  arbitrarily 

choosing  a  vertex  v  ,  partitioning  the  vertex  set  V  into  subsets  V.,...,V  , 

o  1  p 

where  V,  is  the  set  of  vertices  at  distance  d  from  v  (distance  being  measured 
a  o 

by  the  number  of  edges  in  a  shortest  path),  and  finally  deleting  all  the 


edges  joining  vertices  of  the  same  subset  V^,  for  d  =  l,...,p.  It  is  not 
hard  to  see  that  the  resulting  subgraph  is  edge-maximal  triangulated, 
and  its  construction  takes  0(|v|  +  Ie|)  steps.  The  details  follow. 


For  a  graph  G  =  (V,E),  and  a  vertex  of  veV^,  let 
N(v):  =  {ueV j (u,v) eE};  and 

N(v)d+i:  =  N(v)  \  (Vd  U  Vdl)  (with  V_x:  =  0)  . 

W.i.o.g.,  we  assume  that  G  is  connected. 

Algorithm  II  (for  finding  an  EM  triangle- free  subgraph) 

0.  Initialization.  Choose  some  v^.eV  and  set  d  «-  0,  *- 

V  *-  V  \{v^_},  «-  0,  D  <-  0,  and  go  to  1. 

1.  Choose  the  next  vertex.  If  all  veVd  are  marked  and  V  =  0,  stop: 

G  is  an  EM  triangle-free  subgraph  of  G.  If  all  veVd  are  marked  but  V  ^  0, 
set  d  —  d  +  1,  Vd+^  *-  0,  and  go  to  1. 

Otherwise  choose  an  unmarked  veVd,  set  v^  «-  v,  and  go  to  2. 

2.  Add  new  vertices  and  edges.  Generate  N(v  s pt  v  -  v  !' 

- ° —  o  d+1*  aci-  vd+l  d+1 

N(vo'd4.r  v  -  v  0  -  D  1J  (vo-S<''o)d+l)  and  60  to  !■ 

Theorem  2.  Algorithm  II  generates  an  EM  triangle-free  subgraph  in 
0(|V|  +  |E|)  time. 

Proof.  For  any  triangle  of  any  graph  and  any  fixed  vertex  v^,  two  of 
the  three  vertices  of  the  triangle  are  equidistant  from  v^;  hence  the  edge 
joining  them  is  not  included  into  D  by  the  Algorithm.  Thus  G[D]  is  triangu 
lated.  Further,  the  only  edges  of  E  not  in  D  are  those  joining  pairs  of 
nodes  equidistant  from  v^  and  having  a  common  neighbor;  but  adding  any  such 
edge  would  create  a  triangle.  Thus  G[D]  is  maximal. 


To  see  the  complexity  of  the  algorithm  notice  that  for  each 

vertex  v  chosen  in  step  1,  it  takes  0(deg(v  ))  steps  to  identify  N(v  )  ,  , 
o  o  o  d+1 

and  to  add  to  D  the  edge  set  (vq,N(v  )d  ) .  Since  everything  else  done  in 

the  algorithm  takes  constant  time,  the  total  number  of  steps  is  obtained  by 

adding  up  0(deg(v))  for  all  veV,  which  yields  0(|v|  +  |e(),  as  claimed. !i 


Consider  the  complete  4-partite  graph  G  =  K  ,  with 

3  j  3  y  3  9  J 

3  vertices  in  each  of  its  4  independent  sets,  numbered  {1,2,3},  [4,5,6], 
{7,8,9},  {10,11,12}. 

We  first  apply  Algorithm  I  to  G  and  construct  the  EM  triangulated 
subgraph  G[F]  shown  in  Fig.  1.  The  labeling  is  illustrated  in  Fig. 2,  where 
the  circled  numbers  are  the  ranks  (in  the  ordering  a)  and  the  square  brackets 
contain  the  lexicographic  labels. 

Algorithm  I. 

Iteration  1.  Since  all  labels  are  0,  we  choose  arbitrarily  v  =  v,„  -  1. 

n  12 

We  set  a  *-  (1)  and  we  append  12  to  the  label  of  vertices  4,..., 12  (i.e.,  each 
of  these  vertices  now  have  label  L(j)=[12],  j=4,...,12). 

Iteration  2.  We  choose  v^  -  4,  set  a  -  (4,1),  F  -  {(4,1)};  and  append 

11  to  the  label  of  vertices  2,3,7, _ ,12  (i.e.,  we  set  L(2)=L(3)  <-  [11],  L(7)= 

...=L(12)  -  [12,11]). 

Iteration  3.  We  choose  v^0  «-  7  (as  7  is  one  of  the  vertices  with 
lexicographically  largest  label),  set  a*-  (7,4,1),  F  *-  [(4,1),  (7,1),  (7,4)}, 
and  append  10  to  the  label  of  the  vertices  2,3,5,6,10,11,12. 

In  the  next  3  iterations  we  choose  v^  -  10,  Vg  «-  11  and  v^  —  12.  At 
the  end  of  iteration  6  we  have  a  =  (12,11,10,7,4,1)  and  F  =  [(4,1),  (7,1), 
(7,4),  (10,7),  (10,4),  (10,1),  (11,7),  (11,4),  (11,1),  (12,7),  (12,4), 

(12,1)}. 

So  far  every  application  of  step  2  has  resulted  in  the  addition  to  F  of 
all  the  edges  joining  v^  to  its  successors.  At  the  next  iteration  this 


situation  changes. 
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Iteration  7.  We  choose  v,  *-  8 ,  since  vertex  8  has  the  lexicographically 

o 

maximal  label  [12,11,9,8,7].  We  set  <j  «-  (8,12,11,10,7,4,1).  The  first  suc¬ 
cessor  of  8  is  12,  so  we  add  (8,12)  to  F;  but  of  the  remaining  successors  of 
8  only  4  and  1  are  adjacent  in  G[F]  to  12,  hence  only  (8,4)  and  (8,1)  are 
added  to  F,  while  the  remaining  two  edges,  (8,11)  and  (8,10),  remain  in  E\F . 
Thus  we  set  F  -  FU{(8,12),  (8,4),  (8,l)i  and  append  6  to  the  labels  of  all 
unnumbered  vertices  adjacent  to  8  in  G. 

Iteration  8.  v$  -  9,  a  -  (9,8,12,11,10,7,4,1),  F  -  F  lJ  {(9,12),  (9,4), 
(9,1)].  The  edges  (9,11),  (9,10)  remain  in  EVF. 

Iterations  9,10,11,12  produce  a  =  (3,2,6,5,9,8,12,11,10,7,4,1)  and 
the  set  of  edges  shown  in  Fig.  1. 

Next  we  apply  Algorithm  II  to  the  graph  G[EF],  shown  in  Fig.  3,  to  find 
an  EM  triangle- free  subgraph.  We  omit  the  isolated  vertices  1  and  12,  which 
can  be  added  to  the  triangle-free  graph  obtained,  (m  general,  if  G[E  F]  is 
disconnected,  we  apply  the  algorithm  to  every  component  that  contains  a 
triangle) . 

Algorithm  II. 

We  arbitrarily  choose  v^  =  2  and  set  V  ^  {2],  V  '.3 , . . .  ,  11 } . 

Iteration  1.  We  choose  vq  =  2  and  mark  (underline)  2 :  Vq  -  \.2_}.  We 
generate  N(2)^  =  14,5,7,8,10,11}  and  set  *-  N(2)^,  V  —  V  f3,6,9},  and 

D  —  { (2,4),  (2,5),  (2,7),  (2,8),  (2,10),  (2,11)}. 

Iteration  2.  vq  •-  4,  *-  [4, 5 ,7 , 8 , 10 , 11 } ,  N(4)2  {3},  V  ^  —  {3}, 

V  -  {6,9},  D  —  D  'J  {(3,4)}. 

Iteration  3.  v  -  5,  V  -  {4,5,7,8,10,11},  N(5)  -  {3},  D  -  D  U  {(3,5)1. 

O  i- 


Iterations  4  and  5  generate  the  sets  N(7>2  =  N (8)^  =  (3,6),  and  add  to 
D  the  edges  (3,7),  (6,7),  (3,8),  (6,8).  Iterations  6  and  7  generate  N(10)2  = 
N(ll)2  =  {3,6,9},  and  add  to  D  the  edges  (3,10),  (6,10),  (9,10),  (3,11),  (6,11), 
(9,11). 

Iteration  8.  At  this  point  =  {4,5,7,8,10,11}.  i.e.  all  veV  are  marked, 
and  V  =  0;  hence  we  stop,  with  the  graph  G[D]  shown  in  Fig.  4. 

The  original  graph  G  =  K3  3  3  3  had  54  edges.  G[F]  has  30  edges  and 
G [D]  has  15.  Thus  the  edge-maximal  subgraph  G [ F’JD]  with  the  TR-formative 
edge  coloring  [F,D]  has  45  edges. ]j 
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